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\ Abstract: We present a non-perturbative study of the Ac/)^ model in a three dimensional 

\ euclidean space, where the two spatial coordinates are non-commutative. Our results are 

obtained from numerical simulations of the lattice model, after its mapping onto a dimen- 
sionally reduced, twisted hermitean matrix model. In this way we first reveal the explicit C ^ 

' phase diagram of the non-commutative A(/>^ lattice model. We observe that the ordered | p-x 

regime splits into a phase of uniform order and a phase of two stripes of opposite sign, 
^ . and more complicated patterns. Next we discuss the behavior of the spatial and temporal 

correlators. From the latter we extract the dispersion relation, which allows us to introduce 
a dimensionful lattice spacing. To extrapolate to zero lattice spacing and infinite volume 
we perform a double scaling limit, which keeps the non-commutativity tensor constant. 
The dispersion relation in the disordered phase stabilizes in this limit, which represents 
a non-perturbative renormalization. In particular this confirms the existence of a striped 
phase in the continuum limit, in accordance with a conjecture by Gubser and Sondhi. The 
extrapolated dispersion relation also exhibits UV/IR mixing as a non-perturbative effect. 
Finally we add some observations about a Nambu-Goldstone mode in the striped phase, 
and about the corresponding model in d = 2. 
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1. Non-commutative field theory 

The idea of introducing non-commutative (NC) coordinate operators — or a spatial un- 
certainty — dates back to private communications involving Heisenberg, Peierls, Pauli 
and Oppenheimer. The coordinates are then given by hermitean operators obeying a 
commutation relation of the form 

[x^,Xu] = iQlMU ■ (1-1) 

The first papers on this idea appeared in 1947 on quantum theory in flat |j^] and curved Q 
NC spaces. However, it was only at the end of the twentieth century that it attracted 
attention on a large scale in particle physics. Meanwhile, the mathematical foundation for 
quantum field theories on a NC space, NC field theories, was worked out; for an overview 
see refs. [^]. First applications using this type of space as a formalism — rather than a 
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possible description of nature — emerged from solid state physics, see for instance refs. Q. 
Nowadays also the possible reality of a NC space is a subject of intensive research. 

It was string theory that finally boosted this field by identifying strings in certain 
low energy limits with NC field theory ||5| . This was the main reason why NC field theory 
became extremely fashionable over the last years, so that more than 1300 papers have been 
produced on it up to now We also take string theory as a motivation for our study, 
but in this work we investigate a NC field theory as such, without working out possible 
connections to other branches of physics. 

A deep qualitative difference from ordinary (commutative) field theory is the occur- 
rence of a non-locality of the range \/||6|| • This property obviously implies conceptual 
problems, but from the optimistic point of view it may just provide the crucial link to 
string theory or to quantum gravity. In fact, there are arguments that a conciliation of 
gravity and quantum theory induces some sort of non-commutativity under quite general 
assumptions 0. This line of thought would naturally include the time into the non- 
commutativity relation (1.1). However, in that setting the problems related to causality 



and unitarity are especially bad Q. According to refs. Q unitarity is on safe grounds for 
the Minkowski signature, but it is not yet clarified if the transition from a Euclidean to 
a Minkowski signature can be justified. For such reasons, much of the literature excludes 
time from the non-commutativity. This will also be our framework, since the use of a 
Euclidean space is vital for our numerical techniques, which will be explained below. 

In the early days people hoped for yet another possible virtue due to the non-locality, 
namely that it would weaken or even remove the UV divergences of the commutative field 
theories Q, and therefore simplify the renormalization. This hope was badly disappointed: 
first it turned out that only part of the UV divergences are removed whereas others remain. 
In particular the UV divergences in the planar diagrams tend to survive the introduction of 
G [|l^]. What makes the situation much worse, however, is that the remaining commutative 
UV divergences do not just disappear, but they are typically converted into some kind of 
IR divergences. So one ends up with a troublesome mixing of divergences at both ends 
of the spectrum, denoted as UV/IR mixing [11|. A simple intuitive picture of this effect 
based on the uncertainty principle is described for instance in ref. [^]. We are going to 
illustrate the structure of such "mixed" divergences with an example in subsection |1.1| . 
Beyond one loop, perturbation theory does not yet provide any systematic machinery to 
handle this type of divergences, which is unknown in commutative field theory.^ Hence 
adopting a fully non-perturbative approach is highly motivated, and this is the goal of the 
work presented here. 

One issue that immediately arises is the question if UV/IR mixing is a technical prob- 
lem of perturbation theory, or if it should rather be considered as a fundamental property 
of many NC field theories. Our previous numerical investigation of 2d NC U{\) gauge 
theory clearly supported the latter point of view [16|, in agreement with theoretical ar- 
guments [|l^. Manifestations of UV/IR mixing can also be observed beyond perturbation 



""^For completeness we should mention that there are also models known where UV/IR mixing occurs, 
but nevertheless renormalizability can be shown to all orders. Examples for this are the NC Wess-Zumino 
model |l^], the photon self-energy in NC QED [Q, and the real, duality-covariant 4d NC A(^* model ]l5| ]. 



theory — as the present work will demonstrate again — hence such mixing effects belong 
to the very nature of the corresponding NC models.^ This implies, for instance, that it is 
not promising to try to avoid such effects by performing some non-standard perturbative 
expansion (as it has been suggested occasionally in the literature). 

In this work we consider the simplest case of an NC plane with a constant non- 
commutativity tensor, i.e. we deal with the non-commutativity of only two coordinates. 
Its extent is characterized by the parameter t?. 



[xj, Xj] = i'&e 



i,j G {1,2}, 



where e is the antisymmetric unit tensor. In addition we have a Euclidean time coordinate 
t. The whole 3d space is lattice discretized with equidistant lattice spacings. This is 
done by the standard procedure for the time coordinate, while we follow the instruction of 
refs. [17| for the NC plane. There we impose the operator equation 
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2vr ^ 
i — X 
a ■ 
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(1.3) 



where a is the spatial lattice constant. Thus the spatial lattice sites are somewhat fuzzy. In 
our formulation the momentum components pj commute, and we have the usual periodicity 
over the Brillouin zone. 
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(1.4) 



Multiplication by a factor exp(— i ^^^j^p^x^) leads to the condition 
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Hence our NC lattice at fixed "d is automatically periodic. 

Let us assume now that we are dealing with a torus of lattice size N x N: then the 
momenta can take the discrete values 
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l,...,N. 



Together with relation (|l.5| ) we identify the non-commutativity parameter as 

^ = -Na'^ . 
vr 



(1.6) 



(1.7) 



Taking the continuum limit of a NC theory requires to keep i? finite. By condition (1.7), 
the continuum limit a — > and the thermodynamic limit N ^ oo are entangled; we will 
take them simultaneously in such a way that the product Na"^ remains constant. Then 
obviously also the physical volume diverges. This kind of limit is denoted as the double 
scaling limit, and the entanglement that it involves is related to the UV/IR mixing: based 



on relation (1.7) we have to take the UV limit and the IR limit simultaneously. 



^ UV/IR mixing effects have also been observed on tlie semi-classical level ]l^ . 
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1.1 The non- commutative A(/)^ model 

NC field theories can be formulated in a form which looks similar to their commutative 
counterpart, if all field multiplications are performed by the star product [f^ 



{x) -k ■ip{x) := (j){x) exp ( ^'%@f,uOt ) , 



(1.8) 



as we recognize if the fields are decomposed into plane waves; a derivation is sketched in 
appendix |A[ Then we can work with ordinary coordinates x^, and the non-commutativity 
is encoded in the star product. 

In the specific case of bilinear terms in the action, the star product is equivalent to the 
ordinary product, as we see from integration by parts and the antisymmetry of the tensor 
0. Equipped with these rules we can write down the action of the NC Xc/)^ model in the 
d-dimensional Euclidean space, 



(fx 



1 TTl A 

-d^(l){x)df,(l){x) + —4>{xf + -^x) -k (t>{x) -k (f){x) -k <f>{x) 



(1.9) 



We see that the strength of the self-interaction A also determines the extent of the effects 
due to the non-commutativity. 

To illustrate now the UV/IR mixing — that we mentioned before — we consider the 
perturbative expansion of the one particle irreducible two-point function. 



r(p) = (0(p)*0(p))iP/ = ^A^rW(p), 



(1.10) 



where cp is the scalar field in momentum space. To the leading order 1 = the action 
is bilinear, hence the star product is not needed and we find for the free field the same 
result as in the commutative case, T^^\p) = + vn?. However, if we move on to F^^^ 
the commutative term splits into a planar contribution, which is not affected by 0, plus a 
non-planar term, which is altered by the non-commutativity, T^'^^p) = + Fnp (p). 



r(i) 
p 



{2-kY (f' + 



(27r)° 



(1.11) 



The planar term confirms that part of the commutative UV divergences persist. It has 
been shown that this behavior holds for the planar terms to all orders ||l^, 11 1. 

The non-planar term can be evaluated for instance with the Pauli-Villars regulariza- 
tion. In d 



4 one obtains |11, 12| 



1 



967r2 



Ae^fr--^log(^)+0(l) 



A2 



1 



1 



(1.12) 



A is the usual momentum cutoff as it appears in the commutative model and in the planar 
term of eq. ( |1.11 ). Eq. ( 1.12| ) shows explicitly that T^pip) is UV finite (with respect to 
A) in the NC model at finite external momentum p. It diverges, however, if we take in 
addition the IR limit p — > 0, or the limit — > 0. 
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A step beyond standard perturbation theory in the discussion of the NC X(p^ model 
was undertaken by Gubser and Sondhi, who performed a self-consistent Hartree-Fock type 



one- loop calculation |2C]. This method would be exact for the 0{N) symmetric cr-model 
in the large limit. Gubser and Sondhi conjectured its (qualitative) applicability also 
at A = 1. In particular they derived a prediction for the qualitative feature of the phase 
diagram. A strongly negative coefficient corresponds to a very low temperature, such 
that some ordered structure is enforced. 



In dimensions d = 3 and 4, with two NC coordinates obeying eq. (1^), ref. |20| 
predicted that 

• at small ■!?, decreasing tv? results in the usual Ising type uniform order, as in the 
commutative model 

• at larger values of the ordering favors the formation of stripes in some direction, 
where (p{x) takes opposite signs, or even the superposition of stripes in different 
directions ("checker board patterns"). 

The phenomenon of a spontaneous stripe formation is also known in solid state physics, 
see for instance ref. |2l|. However, such a phase is unknown in the commutative X(f)* model. 
Its predicted appearance can be understood as a consequence of an IR divergence, which 
shifts the energy minimum to non-vanishing momenta. At <C the modes close to 
the energy minimum may condense, which corresponds to the stabilization of such stripe 
patterns. 

The conjecture by Gubser and Sondhi was later supported by a study using an effective 
action, which was treated approximately by the Raileigh-Ritz method ||2^. The possible 
existence of this striped phase was also discussed in the framework of a renormalization 
group analysis in 4 — e dimensions . For generalities on the wilsonian renormalization 



group in NC field theory, see refs. [24| 



However, this new phase remained on the level of conjectures until the system could 
be studied non-perturbatively by means of Monte Carlo simulations. The explicit result 
for the phase diagram and for the dispersion relation on the 3d lattice, including its limit 
to the continuum and to infinite volume, will be presented here. Some points of this work 



have been anticipated in various proceeding contributions |25] and in a Ph.D. thesis [26|. 

We note that the occurrence of a striped phase implies the spontaneous breaking of 
translation and rotation invariance. The properties of the Nambu-Goldstone boson related 
to the broken translation symmetry will be discussed in section |^. Based on this symmetry 
breaking, Gubser and Sondhi did not expect a striped phase in d = 2. However, it was 



observed numerically that a striped phase does exist in two dimensions |27], and we are 



going to confirm that observation in appendix As ref. [27| pointed out, this does not 



contradict the Mermin- Wagner Theorem |28|, since the proof for that Theorem assumes 



properties like locality and a regular IR behavior, which do not hold in most NC field 
theories. 



2. The matrix model formulation 



A lattice formulation of the model (|1.9| ) is obtained from the discretization as described 
in section ||, where the derivatives can be replaced (for instance) by differences between 
nearest neighbor lattice sites. However, in such a formulation the field variables on any 
two lattice sites are coupled through the star product. This property would make a direct 
simulation very tedious. 

A way out of this problem is the mapping of the lattice model onto a dimensionally 
reduced, twisted matrix model. This procedure was suggested in refs. as a refinement 



of a previous work on NC gauge theories in the continuum [29|. Some aspects of this map 
are summarized in appendix 

For our scalar field 4'{x, t) on a periodic lattice of size N"^ x T and unit lattice spacing, 
this mapping leads to the action 



s[4>]=nttY, 



t=l 



\ J2 {r,mt] - my + \ {kt + 1) - m)" + "^ktf + \m' 

(2.1) 

where (j){t) represents a hermitean N x N matrix, living on one of the discrete time points 
t = l,...,T. 

Since the time direction is commutative, its contribution to the kinetic term can be 
discretized in the usual way. On the other hand, in the spatial directions the shift by 
one lattice unit has to be arranged for by some matrix transformation, in our case by 
the multiplication with the twist eaters Tj. The condition for them is that they obey the 
't Hooft-Weyl algebra 

r^Fj = ZjiVjVi , (2-2) 

where the phase factor Zji = Z^- is the twist. 

In general a phase factor like Z12 in a twisted matrix model can take the form pO| 

Z12 



^2nik/N ^ ^ g ^ ^ ^2.3) 



since the twist originates from the boundary conditions before compactification. In our 
formulation we set 

*=^, (2.4) 



as we are going to explain in subsection 2.1. Thus our twist differs from the conventional 
choice k = \. Of course, this means that we have to use odd values for N . With this 
choice, the 3d lattice model and the Id twisted matrix model (where the twist appears 
implicitly in the form of the twist eaters) can be rigorously identified [^] by means of 
Morita equivalence [^], which means that the algebras in both formulations are identical. 

Note that typical matrices (j){t) are densely filled. This is in contrast to ordinary 
lattice field theory simulations, where the matrices that appear upon integrating out some 
of the fields tend to be sparse. Hence the matrix multiplications in action (|2.1| ) require 
a considerable computational effort (it grows as 0{N^), i.e. faster than the system size 
of 0{N'^)), since the standard lattice techniques for matrix operations cannot be applied. 



For a description of our numerical methods we refer to appendix j^. At this point we 
only remark that the property of densely filled matrices can be considered as a problem 
inherited from the star products in the lattice action, though the mapping onto the matrix 
model still simplifies the simulation drastically. 



2.1 The twist 

Let us briefly discuss in this subsection our choice of the twist factor. It is instructive to 
return to a general spatial lattice constant a in this context. First, since we are on a torus 
we should not use the unbounded operators xj. Instead we introduce the unitary matrices 



which obey the commutation relation 

4^2 



27r 



(2.5) 



ZiZj 



exp 



iQij ) ZjZi = exp ( -j^i^^ij I ZjZi . 



47r 



(2.6) 



In the last step we have used the relations ( |1.2D and (|1.7|). 

In order to describe the shift by one lattice unit we need the translation operator 



(2.7) 



the operator dj is introduced in appendix To perform this shift correctly, Dj has to 
fulfill the relation 



(2.8) 



The issue is now to find a matrix solution for the conditions ( |2.6D and ( p.Sj ). To this 
end we use the unitary twist eaters 



/O 1 

1 



Fi = 



Vi 



1 
0/ 



/I 



-21 



^21 



^21 



V 



(2.9) 



where Z12 = characterizes the twist. Now we see that the ansatz 



Z^ 



\(N+l)/2 



Zo 



(iV+l)/2 

1 



(2.10) 



does indeed solve the conditions ( |2.6D and ( |2.8D for odd A^, iff the twist has the form 
anticipated in eqs. ( |2.3D and (|2.4|) . Hence this form represents a satisfactory solution to 
the two matrix conditions given above (this solution is unique up to symmetries of the 
algebra) . 



3. A momentum dependent order parameter 



We denote by (j){p, t) the scalar field with the two spatial components expressed in momen- 
tum space (the Fomier transform can be carried out in the usual way, see appendix |A|). 
Now we introduce the quantity 



T 



t=l 



(3.1) 



The field </>(p, t) is averaged over the time direction and then rotated such that its absolute 
value is maximal at a specific value of The expectation value {M{k)) is our momentum 
dependent order parameter. 

Of course, (M(0)) is the standard order parameter for the Z2 symmetry (magnetiza- 
tion). For /c > this order parameter is sensitive to some kind of staggered order, as it 
can occur if anti-ferromagnetic couplings over some distance are involved. In particular, 
(M(l)) detects the formation of exactly two stripes parallel to one of the axes; more pre- 
cisely it measures the leading sine component of ^ 0(x, t) in such a pattern. The value 
k = ^/2 captures the corresponding case with diagonal stripes, and {M{k > 2)) is suitable 
for the search for orders with higher modes of ^ Ylt 'Pi^i^)- 

As an example, figure g shows the results for (M(0)) and (M(l)) at iV = T = 35 
and small A (on the left) resp. large A (on the right). As m? decreases, the uniform resp. 
staggered order is observed to set in unambiguously. 

However, from such plots the critical value of m? can only be read off approximately. 
For a more precise localization of the disorder-order phase transition we measured the 
connected two-point function (M(fe)^)c, which has a peak at the phase transition (for 
— > 00 it diverges at this point). In figure ^ we show the results for the phase transitions 
illustrated before in figure g. Indeed we find marked peaks which allow for an accurate 
determination of the transition from the disordered phase to the uniformly ordered phase 
(on the left) or — at larger A — to a striped phase (on the right). 




Figure 1: The momentum dependent order parameter (M(fc)) (at A; = 0, 1) against N'^m'^ at 
N = 35. On the left we fixed A^^A = 70, which leads to the uniform phase. On the right we set 
N^X = 350, so that a strongly negative leads to the striped phase. 
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Figure 2: The connected two-point function of the momentum dependent order parameter, 
{M{kf)c against iV^m^ at iV = 35 and iV^A = 70 (on the left) resp. TV^A = 350 (on the right). We 
observe marked peaks, which aUow us to locahze weh the phase transitions that we saw before in 
figure |l|. 
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Figure 3: The order parameter {M{k)) at different momenta k, at N = T — 45. 



In figure ^ we take a look at {M{k)) at the (discrete) values of k in the range k = . . .7. 
As long as is not strongly negative, no order can be found at any frequency, so we are 
manifestly in the disordered phase (plots on top). Moving on to <C we find at small 
A clearly the uniform order, since only (M(0)) deviates from zero (plot below on the left). 
Finally at larger A the Ising order parameter (M(0)) drops to zero again, and we observe a 
clear signal at A; = 1, i.e. a two-stripe pattern (plot below on the right). In this case, we also 
see some signals at /c = 3 and 5. However, this should not be interpreted as an indication 
of an underlying multi-stripe structure. What happens is that only close to the phase 
transition the two stripes approximate a sine shape; for even lower m? they approximate 
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Figure 4: The profile $(a;i) = X^ajj t 't'i^) '^^ ^ simple two-stripe pattern at = 35, iV^A = 350: 
near the disordered phase a sine shape dominates (on the left, N'^m? = —250), but at even lower 
m? the profile approximates a step shape (on the right, N'^m? — —620). (These two plots are based 
on single configurations, and the error bars refer to the sum in $(a;i).) 

more and more a step shape, and what we see here is simply the Fom'ier decomposition of 
such a pattern. This transition is illustrated in figure ^. 

A really stable multi-stripe pattern is hard to find — as we will discuss later — and 
in that case we should see {M{k < V2)) ~ and find a signal only at higher k. 

4. The phase diagram 

In the previous section we explained our tools to explore the phase diagram. In particular, 
decreasing at fixed A we could localize the disorder-order phase transition accurately, 
and determine the type of order that emerged. 

The next step is to repeat this procedure at various values of N = T and search for 
suitable axes so that the phase diagram stabilizes for increasing A^. The result is shown in 
figure ^. We obtained an explicit phase diagram with the qualitative features conjectured 
by Gubser and Sondhi. 

It turns out that the axes N'^m?' and A^^A are suitable for a large extrapolation. The 
power of A^ multiplying the self-coupling is singled out by the transition inside the ordered 
regime between the uniform and striped phase. That phase transition, however, cannot 
be localized to the same precision as the disorder-order transition. The reason for this 
property will be clarified later on in the discussion of the dispersion relation (section |6|). 
Still the stabilization of this region in A^ is ultimately compelling. 

Of course, the question about the order of the phase transitions is of interest, but (as in 
many other models) it is difficult to arrive at an absolutely safe answer. To get insight into 
this question, we searched for a hysteresis behavior by crossing the phase transition lines. 
If we cross the transition between disorder and order, we do not see any hysteresis effect. 
We have checked this behavior to a high precision for the transition from the disordered 
phase to both, the uniform and the striped phase, see figure ^. Hence we assume that 
transition to be of second order. 

The situation with respect to the uniform-striped transition is unclear: for increasing 
A in the ordered regime, stripes do eventually show up. On the other hand, for decreasing 
A it seems to be hardly possible to make the stripes disappear again. This behavior is 
illustrated in figure ^. 
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Figure 5: The phase diagram of the 3d Xcj)* model with two NC coordinates, identified from 
(indirect) lattice simulations with lattice size N^. 
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Figure 6: An illustration of the hysteresis behavior described in the text: for the disorder-order 
phase transition, no hysteresis can be observed, as we show for the case that the uniform resp. 
striped order is involved (above, on the left resp. on the right). Below we show the outcome for 
the uniform-striped transformation: once a stripe pattern is built, a return to the uniform order is 
hardly possible inside the ordered regime, hence the "thermal cycle" does not close. 
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Figure 7: Typical lattice configurations at = T = 45 for low — (on top) and high —m^ 
(below), and for weak A (on the left) and strong A (on the right). The dotted and blank area 
represent the domains with different signs of (/'(x). We recognize disorder in both cases on top, and 
a uniform resp. striped order below. 



It is instructive to visualize typical configurations in the different phases. This can 
be achieved by mapping back the matrix configurations to the lattice. The mapping and 
inverse mapping between lattice configurations and matrices are briefly described in ap- 
pendix For further details we refer to ref. ||2^. Typical snapshots at low and at high 
values of A and of —m'^ are shown in figure 0. 

In particular we clearly see a two stripe pattern in the plot on bottom on the right. 
The search for such obvious pictures also for multi-stripe and checkerboard patterns — 
which are safely stable in the Monte Carlo history — turned out be very difficult. Quite 
clear hints for such patterns can be seen at relatively large and a strong self-coupling A; 
examples are shown in figure ^. 

However, if we just map the configurations back to the lattice and illustrate it as in 
figures 1^ and g the emerging pictures do not show a convincing dominance of one stable 
multi-stripe pattern, unlike the case of two stripes parallel to one axis at lower A. 

If we really want to confirm the Gubser-Sondhi conjecture we need evidence that in 
the limit, which corresponds to the continuum and to infinite volume, a finite stripe width 
dominates. Hence in this limit we expect an infinite number of stripes in various directions, 



(a) iV = 35, A = 10, = 
-4 



(b) iV = 55, A = 50, 
= -22 



Figure 8: Examples for a pattern with two diagonal stripes (on the left) and of four stripes 
parallel to an axis (on the right). They were observed in the histories at = 35, A = 10, = —4, 
resp. N = 55, A = 50, = —22. However, it is difficult to verify the ultimate stability of such 
non-minimal stripe patterns. 



with a finite width. Since it is very difficult to find evidence for this behavior by a direct 
illustration, we will search for it in a more subtle manner in sections I and 0. That 
investigation is based on correlation functions, so this is what we want to discuss next. 



5. Correlation functions 

We first take a look at the spatial correlation function^ 

C{x) = ^ + y, t)</'(y, t)) . (5.1) 

As an example we plot in figure ^ the correlators C{xi,Q) and C(0,X2) at = T = 45 in 
the four sectors of the phase diagram, which we also distinguished in figure 1^. In all the 
sectors we confirm the expected behavior: a fast decay in both directions in the disordered 
phase, but hardly any decay in the uniform phase. In the case of two stripes (parallel to an 
axis), we observe a ferromagnetic behavior parallel to the stripes and an anti-ferromagnetic 
behavior vertical to them. 

Next we turn our interest to the fast decay in the disordered phase, at some point close 



to the ordering phase transition. Figure 10 shows logarithmic plots for these correlators, 
one close to the uniform phase (on the left) and the other one close to the striped phase (on 
the right). We see that the decay is somehow irregular: it is faster than polynomial, but it 
does not follow an exponential either. Of course the exponential decay is standard in the 
commutative world, so here we see that this property is distorted by the non-commutativity. 

Now we turn our attention to another type of correlation function, namely to the 
quantity 



G{t) := i Y.^4>{P = 0' ty~^(P = 0. t + (5-2) 



^Theoretically one could also fix y and t, but the summation is useful in practice to enhance the statistics. 





(a) N^X = 90, N^m^ = -22.5 



(h) N'^\ = 900, iV^m^ = -360 





(c) N^X = 90, iV^m^ = -225 



(^rf; iV^A = 900, iV^m^ 



-945 



Figure 9: The spatial correlation function C(x), defined in eq. ( pTl| ) (with one component of x 
being 0) in the different sectors of the phase diagram in figure ||, at = T = 45. In the disordered 
phase (the two plots above) we see a fast decay in both directions, both, for small and for large A. 
In the uniform phase (below, left) there is hardly any decay, and in the case of two stripes (below, 
right) we observe a mixed ferromagnetic and anti-ferromagnetic behavior in the two directions. 
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(a) N^X = 70, N^m^ 



-17.5 



(b) N'^X = 350, iV^m^ 



-140 



Figure 10: The spatial correlator in the disordered phase, close to the transition to the uniform 
phase (on the left) resp. to the striped phase (on the right). The logarithmic plots show that the 
decay is fast, but it does not follow an exponential. 
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This function measures the correlation between two averaged spatial layers with a temporal 
separation of r. We denote G(r) as the temporal correlator, but we should stress that the 
difference from the spatial correlator is not only that we deal with a separation in time. 

In figure |l^ we show typical examples also for this correlator in the four sectors of 
the phase diagram, in analogy to figure Here we show directly logarithmic plots, which 
illustrate that this decay does follow an exponential — resp. a cosh function, i.e. an expo- 
nential with periodic boundary conditions — in the disordered phase and in the striped 
phase. In the uniform phase we obtain G{t) ~ 1 for all values of r. 
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(a) iV^A = 70, N'^m^ = -35 




10 15 20 25 30 35 



(h) N'^\ = 525, iV^m^ = -140 




5 10 15 20 25 30 35 



(c) N'^X^ 70, iV^m^ = -140 




(d) N^X = 525, iV^m^ = -560 



Figure 11: The temporal correlator G{t), defined in eq. (5^), in the different sectors of the phase 
diagram in figure ||, at = T = 35. We see a fast decay, which now follows well an exponential 
(resp. a cosh shape) in the disordered phase (two plots above) and in the striped phase (below, 
right), whereas G{t) is nearly constant in the uniform phase (below, left). 

Prom the exponential decay in the disordered phase close to ordering we can now 
extract the energy at momentum p = 0, i.e. the rest energy Eq. We evaluate this quantity 
at some time separation r as 

G(r + 1) 



En 



In- 



G(r) 



(5.3) 



Such a term is sensible if this evaluation yields a plateau over some range of r, which is not 
very close to the boundary nor to the center T/2, so that the exponential decay dominates. 



Figure 12 shows an example of the corresponding decay and for the neat plateau that one 



obtains for Eq according to eq. (p^)- 

The same procedure to extract the effective energy from an exponential decay can be 
repeated also for finite momenta p. This leads to the full dispersion relation, which is the 
subject of the next section ^. For more general aspects of dispersion relations in NC field 



o 

o 
o 

o 
^^^^ 



theory we refer to ref. |32]. 



6. The dispersion relation 



We now generalize the correlation function G(r) of eq. (5.2) and study the decay of the 
expectation value 

i5](0(p,t)*0(p,t + r)) (6.1) 
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Figure 12: On the left: the temporal correlator G(t) in the disordered phase, close to the transition 
to the striped phase (at iV = T = 35, A^^A = 350, iV^m^ = -140). For comparison, the hne 
represents a cosh function. From the decay we extract the rest energy Eq at different values of r, 



according to eq. (5.3). The plot on the right shows that this provides a stable plateau, up to the 
region around T/2 where the finite T effects are significant. 



at various momenta p. In the disordered phase close to the ordering transition we always 



find an exponential decay in r, which allows us to extract the energy E{p). In figure 13 we 
show the resulting dispersion relation at = T = 35 for a small, a moderate and a large 
value of A. In the first case — near the uniform phase — E'^{p'^) follows closely the usual 
linear form, which confirms again that at small A this model looks like its commutative 
counterpart. 

As A increases we see an amplified rest energy (at p = 0), followed by a sharp dip 
at low momenta and again the linear curve asymptotically at larger momenta. Of course, 
the observation Eq > agrees with the UV/IR mixing. Our result confirms that this 
mixing also occurs non-perturbatively, as we announced in section ^. The question of an 
IR divergence, i.e. a divergence of Eq in the double scaling limit of zero lattice spacing and 
infinite volume, will be addressed in the next section. 

The fact that the energy minimum drives to finite momenta is the basis for the striped 
phase that we illustrated before in section |^. If such momenta condense they manifest 
themselves in the stabilization of some stripe patterns. For instance, in the second plot 
of figure |l^ (at A = 0.6) the leading non-zero momentum which occurs in this volume is 
clearly the location of the energy minimum. This momentum corresponds to A; = 1 in the 
notation of eq. ( |3.lD , hence in this case a condensation of the mode with minimal energy 
leads to two stripes parallel to one of the axes. 



In the last example in figure 13 there are several non-zero momenta which are close 
to the minimum. In such situations various types of stripe patterns may condense. The 
transition between them does not come about easily, so they all appear stable even if they 
do not correspond to the exact minimum. In such cases we could see multi-stripe patterns 
as in figure |8|, but it was difficult to figure out which pattern is ultimately most stable. 
The dispersion relation now clarifies the situation as a coexistence of qualitatively different, 
practically stable patterns in the vicinity of the energy minimum. If we start from such 
a point in the disordered phase and lower m?, the configuration takes one of the stripe 
patterns which correspond to the k values near the minimum. Once this pattern is built 
up, it is hardly possible to change it again as long as we are in the ordered regime (obviously 
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Figure 13: The dispersion relations measured at iV = 35 in the disordered phase, close to the 
ordering transition. We proceed from small A (on top: A = 0.06, w? = —0.01) to moderate A (on 
the left: A = 0.6, — —0.23) and finally to large A (on the right: A — 100, = —7). In the 
first case, which is close to the uniform phase, we see the standard linear dispersion relation. As A 
increases we are close to the striped phase, and the energy minimum moves to a finite momentum 
\p\ . The location of the minimum and the rest energy grow for raising values of A. 
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this would require a large deformation in the structure of the configuration) . For the same 
reason the transition between the uniform and the striped phase could not be determined 
accurately based on a direct illustration, as we described in section ^. However, it is the 
very nature of the system that such transitions happen gradually over a broad interval in A. 

As a further example, we show a dispersion relation for = T = 55 in figure |l^. Here 
we have a finer resolution of the momenta. Now the momentum that corresponds to = 1 
is clearly not the minimum any more, but the minimal vicinity involves the cases k = V2, 
2 and \/5- From the phase diagram in figure H it is clear that for the given parameters 

= — 15, A = 50 we are near the striped phase, so here we have at last an evidence for 
the existence of stable multi-stripe patterns. 



The dispersion relation observed in figure 14 can be parameterized to a good accu- 
racy as 



Co p ^ + + 



ci 



: exp I —C2\Jv'^ + w 



with Co = 0.919(4) , ci = 0.07(1) , C2 = 13(1) , = 0.0004(1) 



(6.2) 



which is illustrated by the curve in figure |T^. (We assume the slight deviation of co from its 
expected value 1 to be a finite artifact.) The square root in this formula is characteristic 
for the three dimensional case, and an obvious ansatz for a linear IR divergence in E'^ 



o 
^^^^ 



-17- 




Figure 14: The dispersion relation determined at iV = 55, ni^ ~ —15, A = 50. The symbols 
correspond to the discrete values of fc = N\p\/ {27r) that occur, i.e. fc = 1, \/2, 2, \/5, VS, ■ ■ ■ 

(regularized by the term fn? and suppressed at larger p^). We remark that in d = 4 
one would expect the corresponding formula with the square root replaced both times by 
(p^ + m^). In d = 3 and 4 this form can be related to Fnp given in eq. (|1.11 ) |p6|| . In our 3d 
results the fits work in fact much better with eq. (|6.2D than with the expected 4d formula. 



7. The continuum limit 

In order to study the extrapolation of our results to a continuous, non-commutative space, 
we first have to identify a dimensionful lattice spacing. We perform this identification in 
the planar limit, where one sends ^ oo at fixed parameters m? and A. 

We saw in section ^ that the dispersion relation E'^{p'^) becomes linear for relatively 



large momenta. In figure |15| we study this behavior in the planar limit: we see that the 
linear regime is not altered any more as we approach this limit, and that the dip for the 
minimum is squeezed towards zero momentum. This implies that the behavior at relatively 
large momenta is dominated by the planar contribution. Hence we can extrapolate a 
dispersion relation to the planar limit by just extending the linear slope down to = 0. 

This is analogous to our study of the 2d NC U(l) gauge theory, where we saw an area 
law for Wilson loops at a small area Thus it coincided with the analytic result by 

Gross and Witten for (commutative) U(A'^) Yang-Mills theory in the planar limit [p^ j. 

Making use of this property, we determine an effective mass Meg in the planar limit. 
We only consider the "linear dispersion regime" where the momentum \p\ is large enough 
to follow the linear dispersion behavior 



M, 



cS+P 



(7.1) 



If the momenta are not too small, this relation is fulfilled to a high precision — as we 



hence we can evaluate an accurate value for in lattice 



saw from figures and |14 
units. In particular, we fixed again A = 50 and considered various values of m^; in each 
case we searched for the planar limit, i.e. a stabilization of the resulting effective mass as 
increases. 
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Figure 15: The dispersion relation (at A = 50, = —13.2) as we approach the planar limit. The 
linear regime is stable, and the minimum moves close to 0. 
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Figure 16: The effective mass squared as a function of the bare mass squared in the disordered 
phase. We observe a linear behavior, which stabilizes accurately in the planar large N limit. 
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Figure 17: The terms fi^ and 7 in the parameterization (7.2), at A = 50. They are found to be 
constant in to a good approximation. 



This is indeed the case in the range = 15 . . . 45. Figure 16 shows that the results for 
Mgg depend linearly on m^, as long as we are in the disordered phase. We parameterize 
this linear curve as 

M^f, = fi^ + jm', (7.2) 

and we show in figure ^ that the values we obtain for fi'^ and 7 are stable in A^. 

Our condition for the continuum limit is now that we keep the effective mass in physical 
units, Mes/a, constant in the planar limit, where a is the lattice spacing. Without loss of 




Figure 18: The spatial correlator C{x) defined in eq. (5.1), amplified by a wave function renor- 
malization factor N^^^^, against the physical distance a\x\. We see a double scaling regime, which 
corresponds to the non-commutative continuum limit. (The dashed line marks the linear regime of 
this double scaling limit.) 

generality we simply set Mefj/a = 1, which means that we measure every quantity which 
has the dimension of a length in units of a/Mcs- This corresponds to identifying the lattice 
spacing as 

(7.3) 



2 2 

H + 7 m 



1^ represents the squared critical bare mass, and according to our 



The ratio 

result illustrated in figure 17 we obtain 



-15.01(8), 



(7.4) 



which is of course also in agreement with figure p^ . 

Therefore, in the double scaling limit — which keeps the non-commutativity parameter 
'& in eq. ( |1.7] ) constant — we have to fix N{fi'^ + 7m^), while taking the limits N ^ oo 
and m^. This limit involves a free constant, which fixes the value of the non- 

commutativity parameter We choose it as 



Na^ = IOO7 



N{m^ - ml) 



vr 



= 100, 
9.77(6) . 



(7.5) 



The value for N{m? 



-mz 



is taken rather large, so that we do not get too close to the critical 



mass for the values that we are using (otherwise we would risk numerical problems due 
to large fluctuations).^ 

In this framework, we first measure the spatial correlator, which was defined before in 
eq. ( ^.1[ ). A double scaling behavior can be observed if the spatial correlator is multiplied 
by a wave function renormalization factor A^", where the suitable power is found to be 
a = 1.35. Then a scaling region shows up, as figure 18 shows. ^ 



*The question if the dimensionless non-commutativity parameter is rational or not is not an issue here, 
because we hand it over to the computer which makes it rational anyhow. 

^The correlator shown here is not normahzed. This is in contrast to sections, where we set C(0) = 1. 
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Figure 19: The dispersion relation at A = 50 and ■& ~ 9.77 (in physical units), with energy and 
momentum expressed in physical units. We see the linear regime starting at moderate momenta, 
a minimum around (p/a)^<0.1, and a rest energy Eq, which increases rapidly with N, see also 
figure |2^. The plot on top provides an overview, and below we zoom the region around the mini- 
mum. This result confirms that the dispersion relation stabilizes in the double scaling limit to the 
continuum and infinite volume, which corresponds to large N in this figure. 

Finally we are equipped to approach the question, which represents the real challenge 
in this context: we now want to verify if the striped phase persists in the continuum 
limit. 

To this end, we consider the double scaling limit of the dispersion relation. The crucial 
question is whether or not the dispersion minimum stabilizes at a finite momentum. If such 
a minimum survives the double scaling limit, it indicates the existence of a striped phase 
in the continuum, where the value of the minimal momentum characterizes the dominating 
stripe width in the infinite, non-commutative space. 

The plots in figures |l^ demonstrate that at finite momentum the dispersion relation 
does indeed stabilize in the double scaling limit. For = 25 . . . 77, a universal minimum is 
found around (p/o)^;^ < 0.1.^ Of course, this number refers to the specific values of "d and 
A that we have chosen. We recall that here we are working in the disordered phase, which 



also holds for figure 19. It would be very difficult to extrapolate to the double scaling limit 



inside the striped phase, because there the finite N boundary conditions have a very strong 
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^^^^ 



®We omit A'^ = 15 because in that case the correlator G(r) decays too fast to extract the energy. Note 
that also A*' = 25 is problematic in this respect. 
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Figure 20: The rest energy in pliysical units as a function of \/N cx 1/a. We clearly recognize a 
linear IR divergence in the continuum limit. 

impact. Hence our strategy is to work in the disordered phase (where the finite effects 
are rather harmless, as figures |9| and |l^ show) close to the striped phase, and extract from 
there information about the striped phase in the continuum, by identifying the modes that 
will condense at lower m^. 

In particular this result shows that the model is non-perturbatively renormalizable. 
In fact the way we identified the lattice spacing amounts to choosing it in such a way 
that the high momentum region of the dispersion relation scales in the double scaling 
limit. What is really non-trivial, therefore, is that the large scaling extends to the low 
momentum region, where it deviates from the planar linear behavior observed at higher 
momenta. This corresponds to the continuum limit of the non-commutative field theory 
at finite 

We remark that in this model also the renormalizability of the perturbative series 
may seem less problematic than in other cases, since the commutative Xcp^ model is super- 
renormalizable. However, to the best of our knowledge the perturbative renormalizability 
of the 3d NC Xcp'^ model has not been shown in the literature. 

At last we take a look at the double scaling of the rest energy Eq . Figure shows that 
it diverges linearly in ^/N oc 1/a. Therefore we do find an IR divergence in the continuum 
limit, in full agreement with the concept of UV/IR mixing — we recall that also the UV 
divergence in this model is linear, as eq. ( |1.11| ) shows. This IR divergence confirms that the 
non-planar terms survive the continuum limit that we take by the double scaling described 
here, and they strongly affect the continuum behavior. 

8. A Nambu-Goldstone mode in the striped phase 

In the striped phase the translational symmetry is spontaneously broken. According to 
the Goldstone Theorem, we expect the emergence of a zero-energy mode in the continuum 
theory. Generalities about Nambu-Goldstone (NG) bosons in NC field theory are discussed 
in refs. Q. 
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Figure 21: The oriented VEV {<^'{x2)) = Ea;i,t('^'(^' *)> against X2 for A = 2 and = 

— 1.04, —1.2, —1.7. Rotations and shifts are carried out before taking the ensemble average, as we 
described in the text. 

On the lattice the NG mode acquires a small energy since the translational symmetry 
is discretized.^ In this section we present numerical results from the striped phase on this 
issue. 

When dealing with the spontaneous breakdown of a symmetry, one has to define the 
vacuum expectation values (VEVs) with some care to make them meaningful. In figure ^ 
(below, right), for instance, we have rotated each configuration before the measurements 
so that the stripe pattern is oriented in one particular direction. More precisely we have 
rotated each configuration such that the momentum p which gives the maximum in M(k) 

— defined in eq. ( |3.1| ) — points to a particular direction. 

The standard VEV {(j){x,t)) vanishes all over the phase diagram. Alternatively we 
now introduce an oriented VEV, which is based on configurations which are all rotated 
and shifted before the measurements. In this section we restrict our attention to the two- 
stripe patterns parallel to an axis. We define the oriented VEV by a rotation so that the 
stripes are vertical to the X2-axis, and in addition by shifts on the rotated configurations, 

(j)'{x, t) = 4>{XI,X2 - Cr, t)rotated , (8.1) 

which maximize the overlap with a characteristic profile. 

For simplicity we use the step profile sign(A^/2 — 3:2), so we determine the shift a by 



maximizing the overlap Y^^^4>'{x,t)sign{N/2 — X2). In figure 21 we show the results for 
the oriented VEV 



($'(x2)) = -^j;(0'(x,t)). (8.2) 



NT 

Xl,t 



The feature is very similar to figure ^ (where we considered single configurations), and it 
demonstrates that the translational symmetry in the X2-direction is broken spontaneously. 



''in the striped phase, the rotational symmetry is spontaneously broken as well. That symmetry is 
also explicitly broken by the boundary condition, unlike the translational symmetry. The breaking by the 
boundary condition may be irrelevant when the wave length of the stripes is much smaller than the spatial 
extent. However, the wave length has to be much larger than the lattice spacing in order to neglect the 
effects of discretization. This makes it very difficult to study the corresponding NG mode by simulations, 
and we do not attempt it in this work. 



Performing a Fourier transformation, we find that Ylti^iP^^)) non-zero for p = 
(0, ^). Deeper in the striped phase we also have condensations for p = (0, ^n2) with 
712 = 3,5,..., which corresponds to the deformation of the leading sine pattern as we 
have seen in figure ^. Note that the pattern is characterized by anti-periodicity in the 
X2-direction with the period N/2, and a constant behavior in the xi-direction. All the 
modes consistent with this periodicity may condense. Due to the condensation of non-zero 
momentum modes, the momentum is no longer conserved in the striped phase. 

In the continuum we can identify the NG mode in the standard manner. We introduce 



a field (f>'{x,t) in analogy to eq. (|8.l|) , i.e. oriented such that the dominant stripes are 
vertical to the X2 axis, and shifted properly. Let us then decompose this field (j)'{x,t) into 
the VEV and the fluctuation as 

(P'{x, t) = {(j)'{x, t)) + (5(/>'(x, t) . (8.3) 

In the case of the assumed stripe pattern, the VEV is not invariant under the shift in the 
X2 direction. The corresponding NG mode is given by the infinitesimal translation in the 
X2 direction as 

6cf>'{x,t)^^yix,t)). (8.4) 

Thus we find that the NG mode has the same Fourier components as the VEV {(j)'{x,t)) 
itself. 

In order to reveal the existence of the NG mode, we repeat the analysis which we 
performed before to study the dispersion relation in the disordered phase. Let us define 
the connected two-point function 



G{p,q;T)c := 



^Y,{^'{p,tr^\q,t+r)) - I^Y.(^\p,tr)\ l ^^(^v,*)) 

t \ t J \ t / 

(8.5) 

Unlike the behavior in the disordered phase, the two-point function can take non-zero values 
also for p ^ q because of the non-conservation of the momentum. From the exponential 
decay with respect to r, we can extract the energy of the intermediate state which couples 
to the operators considered, as we have done in section |6| in the disordered phase. When 
the two-stripe pattern is formed, the subtraction of the disconnected part is needed for 
p = (0, "^n^^) and q = (0, ^n'^^), where n'^^ and n'^^ are odd integers. It is precisely this 
case where the quasi NG mode is expected to appear. 



In figure 22 we plot the energies extracted from the two-point function G{p,p;t)c 
for some low momentum modes against the bare mass squared. In order to see the 
dependence, we compare the results for N = 21 and = 25. 

The energy extracted from the (^,0) mode and the (0,^) mode coincide in the 
disordered phase, but they differ in the striped phase because the rotational invariance is 
broken. The energy of the (^,0) mode and the (^, ^) mode have a minimum near the 
critical point. The minimum values decrease for larger A^, and they seem to approach zero, 
having a cusp-like dip. Thus these modes are expected to have finite dimensionful energy 
in the continuum limit. 



0.4 



0.3 



E 



E 



0.2 



0.1 



0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



"T" 



"T" 



N = 21, n = (0,1) 
TV = 21, n = (1,0) 
AT = 25, i5 = 0,1 
AT = 25, S = (1,0) 



"1 n 



I □ 



-ffl 



J I I L_ 



-r 



N = 21, n = (0,0) 
N = 21, jl = (1,1) V- -i 

TV = 25, ii = (0,0) 1-^^ 

TV = 25, n = (1,1) 1--0-I 



-1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 



-1.7 -1.6 -1.5 -1.4 -1.3 -1.2 -1.1 -1 -0.9 -0.8 -0.7 

^2 



Figure 22: The energies extracted from the connected two-point function, against the bare mass 
squared for TV = 21 and = 25 at A = 2. We show above the results for the modes (0, ^) and 
(2^,0), and below the results for the modes (0,0) and (^, 2^) (with p = 2TTn/N). 

The energy of the (0, ^) mode is small in the striped phase, and it becomes even 
smaller as N increases. At larger N the energy seems to approach zero not only at the 
critical point, but all the way in the striped phase. This behavior suggests the existence 
of the quasi NG mode in the striped phase, whose tiny non-zero energy is due to the dis- 
cretization of the translational symmetry. Our results also support that the discretization 
effects are governed by the ratio of the lattice spacing to the wave length of the stripe 
pattern, which goes to zero at large A'^ in the striped phase. 

The energy extracted from the zero mode (0, 0) does not have a minimum near the 
critical point and it does not decrease as N increases. This is consistent with our obser- 
vation in the disordered phase that the dispersion relation has an infrared singularity, and 
the (0, 0) mode, in particular, has as energy which diverges linearly with the cutoff scale 



as shown in figure 20 



9. Conclusions 

We have presented a comprehensive non-perturbative study of the A(/>^ model in a three 
dimensional space, with two non-commutative spatial coordinates and a commutative Eu- 
clidean time. The non-commutativity tensor is constant. This system is lattice discretized 
and then mapped onto a hermitean matrix model. This mapping is of great advantage for 
our numerical simulations, which would be extremely tedious to run directly on the lattice 
because of the star product. 



We introduced a set of order parameters which are able to detect the uniform or- 
der as well as various types of stripe patterns. Measuring numerically the expectation 
values of these order parameters, we could explore the explicit phase diagram of the lat- 
tice discretized model, see figure |5[ Its structure is consistent with an earlier qualitative 
conjecture by Gubser and Sondhi: at strongly negative some order emerges, which 
is uniform at small A (as in the commutative case), but which is dominated by stripe 
patterns at larger A. The latter is unknown in the commutative Xcp'^ model; we denoted 
it as the striped phase. The phase transition between the disordered phase and the or- 
dered regime seems to be of second order, both, for the uniform order and for the striped 
order. 

We then studied spatial correlation functions, which decay in an irregular way due 
to the non-commutative spatial geometry. On the other hand, the temporal correlation 
functions decay exponentially. That property allowed us to evaluate the dispersion relation 
of the scalar field, in particular in the disordered phase close to the ordering transition. In 
the vicinity of the uniform phase we find the usual linear dispersion, whereas the vicinity 
of the striped phase corresponds to the case of an energy minimum at a finite momentum. 
In fact, the stripe formation means a condensation of such modes of lower energy than the 
zero momentum mode. 

If the momentum is not that small, we enter — also in the vicinity of stripes — into 
a linear regime of the dispersion relation. From the extrapolation of this linear behavior 
we obtained an effective mass, which allowed us to identify a "physical" (i.e. dimensionful) 
lattice spacing a (in the planar limit). This provides a prescription how to take the limit 
to the NC X<l)'^ model in the continuum: we take the double scaling limit, which keeps 
the product Na^ constant, where N is the size of the lattice resp. the matrices. Thus 
we fix the non-commutativity parameter -d. Approaching this limit, which is at the same 
time an infinite volume limit, we found stable dispersion relations. In particular, the 
shift of the energy minimum to a finite momentum survives this limit, which shows that 
the striped phase is indeed there to stay in the continuum. This limit completes the 
confirmation of the conjecture by Gubser and Sondhi. Moreover, the convergence of the 
dispersion in the double scaling limit shows that this model can be renormalized non- 
pert urbatively. 

The stripe formation implies the spontaneous breaking of translation and rotation 
invariance in the spatial plane. The Nambu-Goldstone mode of the broken translation 
symmetry can be seen approximately from the numerical data. Interestingly, this symmetry 
breaking also occurs in the two dimensional model, see ref. and appendix |^. This is 
possible because this model does not obey the assumptions for the proof of the Mermin- 
Wagner Theorem. 

The energy of the zero mode itself (rest energy) diverges in the limit of zero lattice spac- 
ing and infinite volume. This confirms that the UV/IR mixing is a fundamental property 
of NC field theories, which is also manifest beyond perturbation theory. Diagrammatically, 
this mixing effect is driven by the non-planar diagrams, which are not suppressed in the 
double scaling limit (in contrast to the planar limit). This mixing is the basis for the 
different types of long ranged orders that occur in this model. 



As an outlook, we are now elaborating a dispersion relation for the photon (or, strictly 
speaking, for the corresponding "glueball") in a four dimensional NC space [^]. We have 
studied earlier QED in a NC plane [^], where we already saw striking manifestations 
of UV/IR mixing. We expect such effects also in four dimensions, in particular a 0- 
deformed dispersion relation for the photon. This dispersion has been studied before to one 
loop in perturbation theory [36|, though with an uncontrolled behavior of higher orders. 



In contrast, explicit non-perturbative results would allow us to establish bounds on the 
magnitude of the non-commutativity in nature. 

For instance, blazars (highly active galactic nuclei) are assumed to emit bursts of 



photons simultaneously, which cover a broad range of energies, see e.g. ref. |37]. Experi- 
mentalists are already trying to detect a relative delay of these photons upon their arrival 
on earth, depending on the frequency ||3^ . Such experimental efforts will be intensified 
in the near future. For instance, the GLAST project is scheduled to be launched in 
2006 and to monitor gamma rays from 20MeV up to 1 TeV. If we arrive at results for the 
0-deformed photon dispersion — analogous to figure ^ — they could then be confronted 
with such experimental data p0[. 
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A. The mapping between the lattice and the matrix formulation 

To describe the connection between the lattice formulation and the matrix model that we 
actually simulated, we first go back to the continuum and introduce the terminology of 
Weyl operators. 

We start from a scalar field (j){x) in the (commutative) euclidean space R'^', which falls 
off fast enough for the Fourier integral 



(t){p) = j <fx(t){x)e-'''^P^^ (A.l) 

to converge. The corresponding NC space is characterized by the coordinate operators 
x^, which obey the non-commutativity relation ( |1.1| ). We recall that we are dealing with 
commutative momenta in the NC space. Now we introduce the Weyl operators 



W[<P] = ^ j d^pm e-Mf'^ , (A.2) 

which build a NC, associative algebra. In coordinate space, the map between the scalar 
fields and the corresponding Weyl operators reads 

W[(j)] = J (fx4){x)k{x) , A(x) = j d'^pei(^^-^^)P^ , (A.3) 

where A(x) is a hermitean operator. 



Derivatives of the Weyl operators are defined by an anti-hermitean operator d^, which 
obeys 

[d^,,x,]=5^, ^ [a^,A(x)] = -a^A(x). (A.4) 
This imphes the desired property 

[d^,Wm=W[d^.ct>\. (A.5) 

Hence a translation can be expressed by unitary operators exp(?;^9^), (c.f. eqs. ( ^^ ) 
and (ID), 

K{x + v) = e^"^"A(x)e-''^^^ . (A.6) 

Accordingly a trace on the algebra of Weyl operators is translation invariant, and it is 
given explicitly by 

Tr #[(/>] = j (fx (t){x) , Tr A(x) = 1 . (A.7) 
If G is invertible — which requires the dimension d of the NC space to be even — it 



can be shown that [12| 

Tr(A(x)A(y)) =<5W(x-y). (A.8) 
Thus we arrive at the inverse map 

(/>(x) = TV {w[^] k{x)^ . (A.9) 
Based on this inverse map, we obtain the product of Weyl operators as 

W[(t)\ W[il)] = W[(^ * V] , (A.IO) 



which justifies the use of the star product in eq. (|l.8| ) and below. 

Let us now focus on d = 2 and reduce the plane to a periodic N x N lattice. Then 
also the Weyl operators are reduced to N x N matrices, for instance 

A{x)= Zi"i4"2gxp(-^(nin2 + n^x^)j , (A.ll) 

"1,12 = 1 

where the operators Zj are defined in eq. ( |2.5| ) (here we set a = 1). 

This provides a map between a lattice field and the Weyl matrices, which is analogous 
to the map from continuum fields to Weyl operators, 

X 

The lattice field in momentum space, 

'^(^) = ^ E , (A.13) 

X 

is related to the Weyl matrix as 

= ^(^(n)J(n) , 4>{n) = ^ Tr (j> j'' {7 
n 

where 

J(n) = Zf^is"" e-2™i"2/iv _ (-^ -^4^ 



B. The non-commutative A(/)^ model in d = 2 



In this appendix we want to address the case of only two NC spatial coordinates, i.e. we 
omit the time direction now (resp. we reduce it to one point). 

In this context, we would like to repeat that the formation of stripes necessarily implies 
the spontaneous breaking of translational and rotational invariance. For this reason, a 
striped phase was originally not expected in d = 2 . 

However, a numerical study by Ambj0rn and Catterall revealed that the striped phase 



does exist even in the 2d model |27|. They also mapped the lattice formulation onto 
a matrix model and simulated = 39 at a fixed value of A, probing one line for m^. 
They saw the two stripe pattern along with more complicated multi-stripe patterns. They 



pointed out that this does not contradict the Mermin- Wagner Theorem |28], because the 
derivations of the latter are all based on assumptions like locality, which are not realized 
in this case. Intuitively one may wonder why the short-ranged non-commutativity can 
induce a long-range order. However, keeping in mind that we found UV/IR mixing to be a 
fundamental property also beyond perturbation theory, this does not seem surprising any 
longer. 

We verified the results reported in ref. |27| for d = 2, and extended also that investiga- 



tion to a full exploration of the phase diagram. First we run the simulation with the same 
algorithm as in ref. ||2^: in a Metropolis step the matrix cj) may be replaced hy 4> + erj, 
where is a hermitean random matrix and e is a small parameter which is tuned for an 
acceptance rate around 50 %. After about 500 steps we saw a rich structure of patterns 
in the striped phase, in agreement with ref. |27|. Some typical multi-stripe patterns that 



appeared in this way are shown in figure 23. However, when we extend the history to 
O(IO^) . . . 0(10"^) steps only the two stripe patterns survive, as in the three dimensional 
case. 

We then checked this result with the more powerful algorithm that we also used in the 
3d case. The difference is that now the matrices (p are updated by changing only pairs of 
conjugate matrix elements (pij = (j)*^. This takes more steps, but the corresponding param- 
eter e can be taken much larger. In total this results in a reduction of the thermalization 
time by a factor ~ 1/A^. For further numerical details we refer to appendix |C|. 
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Figure 23: Three snapshots of meta-stable multi-stripe patterns in the 2d model, at iV 
jy3/2j^2 ^ _2i^ and TV^A = 350, in the striped phase. 
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Figure 24: The phase diagram of the 2d NC X(f>'^ model, identified on a lattice of size iV^ . Note 
that the vertical axis has to be scaled differently from the 3d case shown in figure 

The result was that we could see now with less steps the behavior described above, 
namely that in this magnitude of N only two stripes parallel to one axis are ultimately 
stable. Of course, also the variety of meta-stable patterns found first by Ambj0rn and 
Catterall is interesting in itself and their observation complements the description of the 
2d model. 

The feature of the phase diagram is similar to the 3d case, but as a qualitative differ- 
ence the mass axis has to be scaled as N^/'^m? (rather than N'^m?). Then we find again 
an accurate transition line between the disordered phase and the ordered regime, and a 
somewhat broad but stable transition region between the uniform phase and the striped 
phase, see figure |4[ 

Finally we would like to add that the 2d A(/)^ model has recently been studied nu- 



merically on a "fuzzy sphere", which also represents a NC manifold. Ref. |41] observed a 
"matrix phase", which is likely to coincide with the striped phase. For another numerical 
study of a matrix model involving the fuzzy sphere, see ref. p^]. 
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C. The simulations 



Our Monte Carlo simulations were performed with the standard Metropolis algorithm. The 
update steps of a configuration run over the single matrix elements (pij{t) {i > j) which 
may be modified as 

<pij{t) ^ (^ij{t) + e-q , 4>ji{t) ^ 4>jiit) + er/* , e G M+ . (C.l) 

The matrix elements 4>ij{t) = 4'ji{t)* of the initial hot configurations, and the variable are 
complex random numbers, where the real and the imaginary part have a flat probability 
density in [—0.5, 0.5]. The parameter e, which controls the step size, is set in the beginning 
to N~'^l^. It is then adapted dynamically as follows: after updating the complete lattice 
configuration, it is increased resp. decreased by 20 % if the Metropolis acceptance rate is 
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N 




15 


25 


35 


45 


55 


number of configurations 


phase diagram 


300 


180 


60 


30 




2-point function 


10000 


10000 


5000 


3000 


2700 



Table 1: Our statistics at the different values of TV = T in the measurements for the phase diagram 
in figure ||, and for the connected 2-point function (M^(fc))c (examples are shown in figure §). 



larger than 0.6, resp. less then 0.3. Thus the acceptance rate is pushed into this interval, 
and the parameter e typically stabilizes around one order of magnitude below its starting 
value. 

We compared this updating procedure to an alternative one which updates the whole 
matrix 4>{t) at once, c.f. appendix ^. We found the method described above to be superior 
by far with respect to the thermalization time and to the autocorrelation. 

For our method, the autocorrelation time ranged from about 10 to 90 update steps, 
depending on the point in the phase diagram. To keep this effect under control, we pro- 
ceeded as follows: for each simulation we performed six completely independent hot starts. 
Then we thermalized over about 500 to 1000 update steps. After this we measured the 
observables on configurations which were separated by 100 steps each time. 

The statistical error was evaluated with the jack-knife and the binning method. In 
both cases we varied the bin size and we took finally the over-all maximum as a careful 
estimate for the error bar, which is shown in the plots. For the bulk of the points in the 
phase diagram we arrived at small errors already with a modest statistics. To compute the 
2-point functions of the order parameter, however, we had to handle in particular the region 
close to the phase transitions, which required a much larger number of measurements. Our 
statistics is given in table |l|. 

As a test for the correctness of our code, we expanded the expectation value of the 
action anal5d;ically to the first order in A. The result is in excellent agreement with our 
numerical data up to A ~ 0.5. For the corresponding plot, and for further details about 
the simulation, we refer to ref. |26|. 
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